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Summary. — A review of recent nucleon structure calculations within lattice QCD 
is presented. The nucleon excited states, the axial charge, the isovector momentum 
fraction and helicity distribution are discussed, assessing the methods applied for 
their study, including approaches to evaluate the disconnected contributions. Re- 
sults on the spin carried by the quarks in the nucleon are also presented. 



1. Introduction 

During the last decade, results from simulations of QCD have emerged that already 
provide essential input for a wide range of strong interaction phenomena as, for example, 
the QCD phase diagram, the structure of hadrons, nuclear forces and weak decays. In 
this presentation, we focus on hadron structure calculations using state-of-the art lattice 
QCD simulations [1, 2, 3]. 

Understanding nucleon structure from first principles is considered a milestone of 
hadronic physics. A rich experimental program has been devoted to its study, starting 
with the measurements of the electromagnetic nucleon form factors initiated more than 
50 years ago. Reproducing these key observables within the lattice QCD formulation is a 
prerequisite for obtaining reliable predictions on observables that explore physics beyond 
the standard model. In particular we discuss three fundamental quantities, which can 
be easily extracted from nucleon matrix element calculations on the lattice, namely the 
nucleon axial charge, the isovector momentum fraction and the helicity moment. The 
reason for this choice is that they are known experimentally and they all can be directly 
determined at momentum transfer squared q 2 = 0. Thus there is no ambiguity associated 
with having to fit the g 2 -dependence of the form factor (FF), as for example, in the case 
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of the anomalous magnetic moment where one needs to fit the small ^-dependence of 
the magnetic FF. In addition, only the connected diagram, for which established lattice 
methods exist for its calculation, contributes. 

For the discussion of the spin content of the nucleon however, one needs to calculate 
the contributions carried by the u- and d- quarks. Contributions from fermion loops need 
to be taken into account in this case. There are on-going efforts to compute these so 
called disconnected contributions, which will be briefly discussed here. However, in the 
final results shown in this work for the quark spin these contributions will be neglected. 



2. Methodology 

It is convenient to decompose the nucleon matrix elements, (N(p', s')|Op A ' 1 "' M ™ \N(p, s)), 
into generalized form factors (GFFs) as follows 



u N (p ,s ) 



(1) +mod(n, 2)C n+lfi {q 2 )±q^q^ . . . q^ 



u N {p, s). 



For p' = p one has the n th moment of the unpolarizcd parton distribution, {x n ) q . A 
similar expression can be written for QWi—Hn j n terms of A n i(q 2 ) and B n i(q 2 ). The 
ordinary nucleon form factors are then special cases of GFFs given by 



Fi(q 2 ) = J^dxH(x,^,q 2 ), B w (q 2 ) = F 2 (q 2 ) = J^dxE(x,^,q 2 ) 



f l z-i 

„2\ _ n f~2\ 



(2) A 10 (q 2 ) = G A (q 2 ) = J ^dxH(x^,q 2 ), B w (q 2 ) = G p {q 2 ) = J ^dxE{x,^,q 2 ) 

while A n o(0), j4 n o(0), A^ Q (0) are moments of parton distributions, e.g. (x) q = ^20 (0) and 
(^)aij = ^2o(0) arc the first moments of the spin independent and hclicity distributions. 
With knowledge of these moments, one can evaluate the quark contributions to the 
nucleon spin using the decomposition J q = ^[A 2 o(0) + B 2 o(0)} = ^AE + L q and AE = 
Aio. The nucleon spin sum rule \ = ^AE + L q + J g allows the determination of the 
gluon contribution J g . 

In lattice QCD. one proceeds with the calculation of nucleoli matrix elements, by 
creating a state with the quantum numbers of the nucleon at an initial time, U = 
(source), inserting a current at some intermediate time, t, and subsequently annihilating 
the nucleon at the final time, tf, (sink). Such a three-point function is shown in Figs. 1 
and 2, where the lines represent fully dressed quark propagators. For a general operator 
Or, one is interested in the momentum-dependent matrix element: 

(3) G ^-^^-t f -t-q) = £ r^(J«(f / ,t / )O^'-^( a -,t)4(0,0)}e-^ 

Xj,X % 

where the projection matrices T v are To = j(1+7o)j = i^olslk, k = 1,2, 3. Depend- 
ing on the quark structure of the operator Or, such a matrix element is a combination 
of connected and disconnected contributions, as shown in Figs. 1 and 2 respectively. 
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Fig. 1. - Connected contribution to a baryon 
3-point function. 




Fig. 2. - Disconnected contribution to a 
baryon 3-point function. 



The nucleoli matrix element is obtained in the limit where the time separation between 
the current insertion and the source and the sink are large enough such that only the 
ground state matrix element dominates. An appropriately defined ratio 



' G(0,t f ) V G(0,t f -t)G(q,t)G(q,t f ) */-*>i 

thus becomes time independent ( "plateau region" ) and unknown overlaps of the nucleon 
state with the trial states are canceled. In Eq. (4) we have taken the final momentum 
j? = 0. The two-point function appearing in Eq. (4) is given by: 

(5) G(q,t f ) =Y,T^(JUSf,t)4(6,0))e- i ^f. 

Energies of hadronic states are extracted from two-point functions by taking the ratio 
E e s(t) = log G(0, t)/G(Q, t + 1) — » E . Utilizing a basis of interpolating fields Jn(x) 

one can extract both the ground state and excited states. In Fig. 3 we show an example 
of utilizing a five-dimensional basis of interpolators to extract the excited states of the 
nucleon in the positive parity channel. As can be seen, identifying a plateau region is 
crucial to reliably extract the energy of excited states. Using more suitable interpolating 
fields can improve the plateaus as shown in Fig. 4 for both the positive and negative 
parity channels. 

Unlike the case of the isovector nucleon matrix elements, where quark disconnected 
contributions cancel, for the isoscalar case, both connected and disconnected quark dia- 
grams contribute. These contributions are shown in Fig. 5, for the case of the isoscalar 
nucleon matrix element of the axial- vector current at q 2 = 0. Results for the ratio of 
Eq. (4) are shown for various source-sink time separations. When the ground state domi- 
nates there is convergence of the results obtained for various source-sink time separations. 

The disconnected diagrams are notoriously difficult to compute because they exhibit 
large statistical errors and require knowledge of all elements of the quark propagator. The 
so called all-to-all propagator is impractical to compute and store explicitly. Stochastic 
techniques are typically employed to obtain estimates of the all-to-all propagator com- 
bined with techniques to efficiently multiply the number of measurements to suppress 
the statistical uncertainty Such calculations have been shown to be especially suited for 
GPUs [4, 5]. In Fig. 6 we show the disconnected contribution to the isoscalar nucleon 
matrix element of the axial- vector current for various source-sink separations. Note the 
much larger errors as compared to those obtained from the connected part. 
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Fig. 3. - The four lowest states of the nucleon 
in the positive parity channel extracted using 
Nf — 2 twisted mass fermions (TMF) with 
pion mass 308 MeV. The solid lines and bands 
show the value and error extracted from a fit 
in the plateau region. 
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Fig. 4. - The two lowest states of the nucleon 
in the positive (lower) and negative (upper) 
parity channels using two dynamical flavors 
(Nf — 2) of TMF gauge ensembles with pion 
mass 308 MeV. 



3. Results 

The spectrum of the nucleon is explored by several groups. In Figs. 7 and 8 we 
show recent results on the first excited states of the nucleon [6]. All lattice results are 
in nice agreement for the ground state. Results for the first excited state show larger 
statistical errors and deviations. Most lattice data are higher than the mass of the Roper, 
even when pion masses at almost the physical value are used. This is currently under 
investigation. In the negative parity channel, lattice results give the experimental value 
of the (1535). 

Reproducing the nucleon axial charge is important for benchmarking our lattice QCD 
formulation beyond reproducing masses of low-lying hadrons. We show in Fig. 9, the ax- 
ial charge of the nucleon gA, extracted from fitting to the plateau region of the ratio of 
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Fig. 5. - The nucleon matrix element of the 
isovector (isov) and the connected contribu- 
tion to the isoscalar (isos) axial-vector opera- 
tor. The solid line shows the fit to a constant. 




Fig. 6. - The disconnected contribution to the 
nucleon matrix element of the isoscalar axial 
operator as a function of the insertion time- 
slice. The solid line shows the fit to a con- 
stant. 
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Fig. 7. - The positive parity states of the nu- 
cleoli using twisted mass and clover fermions 
from various groups. The magenta filled cir- 
cles show the experimental values of the nu- 
cleon mass and the Roper. 



1.UI 1 1 1 1 1 1 r= — 

■ Twisted Mass (this work) BGR J = A 

□ Clover (this work) # Experiment: N"(1535) 

3.5- ♦ CSSM O S-wave:n+N 




t 



05 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 
ml [GeV] 

Fig. 8. - The negative parity states of the nu- 
cleon using twisted mass and clover fermions 
from various groups. The filled magenta circle 
shows the 5n(1535), whereas the open ma- 
genta circle the ttN mass. 



Eq. (4) and which uses the isovector current for which disconnected contributions can- 
cel. The first observation is the consistency between different lattice formulations. This 
consistency is non-trivial, since the results shown are at different, non-zero lattice spac- 
ings, and in general these different formulations have different finite-spacing correction 
terms. Furthermore, all results underestimate the experimental value as the pion mass 
approaches its physical value. Candidates for this discrepancy are currently thought to 
be systematic effects in lattice calculations, such as artifacts caused by the finite volume 
or contamination of the matrix element by excited states. The latter has been investi- 
gated thoroughly in a high statistics calculation [7]. As can be seen in Fig. 10, qa does 
not suffer from excited state contamination since the plateau value remains consistent 
for a large range of source-sink separations. This conclusion is consistent with Fig. 5. 
An alternative to fitting the plateau of Eq. (4) is the summation method that modifies 
the excited state contributions [8, 9]. Further sources of systematic uncertainties are 
currently under investigation by various lattice QCD collaborations. 
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Fig. 9. - Lattice results on the axial charge of 
the nucleon for: Nf=2 and N f =2 + 1 + 1 
TMF [10]; Nf = 2 + 1 DWF [11]; hybrid [12]; 
N f =2 clover [13]; Nf=2 + 1 clover [14]. The 
physical point is shown by the asterisk. 
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Fig. 10. - The value of gA versus the sink- 
source time separation (Nf — 2 + 1+ TMF at 
~350 MeV. 
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Fig. 11. — Lattice results on the isovector mo- 
mentum fraction {x} u -d of the nucleon. The 
crosses, crossed and open squares are from 
Refs. [15, 16, 17], respectively. The rest of 
the notation is the same as that in Fig. 9. 



Fig. 12. - Lattice results on the isovector po- 
larized moment (x)Au-Aci of the nucleon. The 
notation is the same as that in Fig. 11. 



In Fig. 11 we show the isovector combination for the momentum fraction ^2o(0) = 
(x) u —d for various fermion discretization schemes. The general conclusion is again that 
lattice results are in agreement among themselves, a strong indication that finite lattice 
spacing corrections are smaller than the statistical errors. Results at almost physical 
pion mass in this case appear to be converging to the physical point. Confirmation of 
this trend is on-going with computations being performed at near physical pion mass 
and large volumes. In Fig. 12 we show results for the helicity ^2o(0) = (x)au-AcI- The 
situation is very similar to the momentum fraction. 

The matrix elements presented so far are directly related to the contribution of the 
nucleon constituent quarks to its spin, as explained in the introduction. Knowledge of 
both isovector and isoscalar combinations are required in order to obtain the individual 
quark contributions. From here on, we ignore disconnected contributions to the isoscalar 
matrix elements. There are indications that their contributions are small, e.g. for the 
case of qa = -<4io(0) these are expected to be of the order of 10%, as can be seen from 
Fig. 6. Calculations are underway for their precise determination. The total angular 
momentum carried by quarks requires knowledge of i?2o(0), which on the lattice cannot 
be obtained at zero momentum transfer. i?2o(0) is therefore calculated by fitting the 
momentum dependence, for which we assume a linear form. 

We show in Fig. 13 the quark contributions to the proton spin. We compare results 
obtained using Nf = 2 and Nf = 2 + 1 + 1 TMF with those using a hybrid action of 
Nj = 2 + 1 staggered sea and domain wall valence quarks [12]. As can be seen, the total 
spin carried by the u- and d-quarks J u+d ~ J u since J d ~ 0. This can be traced to the 
intrinsic spin ^AE d and angular momentum L d of the d-quarks canceling. 

To obtain some insight on the spin carried by the quarks at the physical point we 
use Heavy Baryon Chiral Perturbation Theory (HBxPT) to extrapolate the lattice data. 
The expressions for these fits are given in Ref. [10], where also covariant %PT was used 
giving qualitatively similar results. A source of systematic error may came from the fact 
that i?2o(0) is obtained via a fit. This is done by extrapolating the data obtained by using 
two fit ranges in the fit to B2o(Q 2 ), shown by the red and green bands in Fig. 14. As can 
be seen both fit ranges yield consistent results. For the case of the intrinsic spin where 
experimental data are available, a comparison shows that the down-quark contribution 
is under estimated, while for the u-quark there is agreement. However, this disagreement 
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Fig. 13. - Quark contributions to the proton spin as a function of the pion mass. Left: the 
total contributions of the up- and down-quarks (J u and J d respectively). Filled squares are 
obtained using Nf = 2 + 1 + 1 TMF, filled circles using Nf = 2 TMF, compared to results 
obtained using a hybrid action [12] (open symbols). Right: The quark spin and quark orbital 
angular momentum contribution to the nucleon spin obtained using Nf = 2 + 1 + 1 twisted mass 
fermions (filled squares and diamonds respectively), using Nf = 2 twisted mass fermions (filled 
circles and triangles respectively) and using a hybrid action (open circles and open triangles 
respectively). Asterisks denote the experimental results. 



is too small to account for the missing 50% of the spin of the nucleon. Disconnected 
contributions to the quark spins are also expected to be small and therefore a large 
fraction of the proton spin is not carried by the quarks. 



4. Conclusions 

A review of recent lattice QCD nucleon structure calculations is presented, comparing 
the results obtained by different groups. The nucleon first excited states are presented 
with good agreement with experiment in the negative parity channel but still not being 
conclusive in the case of the Roper. Lattice QCD data on the nucleon axial charge are 
in agreement among themselves but, however, underestimate the experimental value for 
qa believed to be due to systematic effects appearing in lattice computations. This calls 
for an in depth study to identify possible sources of this discrepancy, something which 
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Fig. 14. - Chiral extrapolation of the quark contributions to the proton spin. The vertical line 
denotes the physical point. The data are obtained using Nf = 2 twisted mass with a = 0.089, 
0.07 and 0.056 fm (red circles, green triangles and magenta diamonds), and Nf = 2 + 1 + 1 
twisted mass with a = 0.086 and 0.066 fm (blue inverted triangles and blue squares). 
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is being investigated currently by several lattice groups. Lattice data for the momentum 
fraction and helicity on the other hand are seen to converge towards the experimental 
value as the pion mass decreases below 250 MeV although still remain higher. Recent 
results near the physical point have suggested excited states [9] to be the source for the 
remaining discrepancy. This is being investigated. 

The axial charge, combined with the momentum fraction, give insight on the fraction 
of the nucleon spin carried by quarks. Using both isovector and isoscalar combinations of 
these quantities, with a cautionary remark on the omission of disconnected diagrams that 
contribute to the latter, we extract the individual up- and down-quark contributions to 
the proton spin. Qualitatively, we find agreement with what is observed experimentally 
for the quark intrinsic spin fractions. Recent lattice calculations also agree in indicat- 
ing that the total spin fraction carried by down-quarks is much smaller and near-zero 
compared to the up-quark and that the quarks account for about half of the proton spin. 
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